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The application ol modern layer-by-layer growth techniques to transition metal oxide materials 
raises the possibility of creating new classes of materials with rationally designed correlated electron 
properties. An important step towards this goal is the demonstration that electronic structure can be 
controlled by atomic composition. In compounds with partially occupied transition metal d-shells, 
one important aspect of the electronic structure is the relative occupancy of different d-orbitals . 
Previous work has established that strain and quantum confinement can be used to influence orbital 
occupancy. In this paper we demonstrate a different modality for orbital control in transition 
metal oxide heterostructures, using density functional band calculations supplemented by a tight 
binding analysis to show that the choice of non-transition-metal counterion X in transition-metal 
oxide heterostructures composed of alternating LaNiOz and LaXOs units strongly affects orbital 
occupancy, changing the magnitude and in some cases the sign of the orbital polarization. 

PACS numbers: 75.70.Cn, 73.20.-r, 71. 15. Mb, 71.10.-w 



I. INTRODUCTION 

Recent progress in growth of transition metal ox- 

I ide multilayers with atomic-scale chemical precisior^^ 
suggests that it may become possible to create new 
classes of materials with desirable electronic proper- 
ties based on aspects of correlated electron physics 
such as high Curie temperature ferromagnetism^i^, 
'colossal' magnetoresistance^, correlation-driven metal- 
insulator transitions^ and high transition temperature 

I superconductivitji^. The ultimate goal is 'materials by 
design', in other words to construct systems with desired 
electronic properties. A necessary first step is to design 
and fabricate systems with a desired electronic structure. 

The 'correlated electron' properties of transition metal 
oxides are controlled in part by the relative occupancy 
of the different transition-metal d orbitalsii Controlling 

] the orbital occupancy by materials design is therefore an 
important milestone in the progress towards a rational 
design of correlated electron materials. A difference in 
relative occupancies of orbital states may be described 
as an 'orbital polarization' in analogy to the difference in 
occupancies of spin states which gives rise to spin polar- 
ization. As the control of spin polarization is achieved by 

[ application of appropriate magnetic fields, so the control 
of orbital polarization may be achieved by identification 
and manipulation of appropriate 'orbital fields'. 

Two classes of 'orbital field' are well established: strain 
and quantum confinement. Lattice strain works because 
in a transition metal oxide the hybridization between the 
transition metal d orbital and the oxygen p states pro- 
duces ligand fields whose strength depends on the geo- 
metrical distance between the transition metal and oxy- 
gen site. An applied strain changes relative bond lengths, 
thereby affecting the ligand fields. Strain control was re- 
ported by Konishi and collaborators^, who grew films of 
'colossal magnetoresistance' (CMR) manganitcs on sub- 



strates which imparted compressive, negligible, or ten- 
sile strain to the manganite film. The systems exhibited 
resistivities which were respectively strongly insulating, 
weakly insulating, and metallic and the change in resis- 
tance was attributed to a strain-induced change in rel- 
ative occupancies of the d^^-y^ and d^^^-r^ orbitals on 
the electrically active Mn site. 

The quantum confinement effect works because each 
transition metal d-orbital has a specific spatial structure 
which leads to a direction-dependent hopping amplitude. 
Spatially anisotropic quantum confinement (for exam- 
ple in a heterostructure composed of alternating tran- 
sition metal oxide and insulating spacer layers) allows 
electrons to delocalize more in some directions than in 
others, thereby allowing some orbitals to gain more de- 
localization energy than others. The quantum confine- 
ment effect operates in a slightly subtle way in transition 
metal oxide systems, where the relevant electronic bands 
are antibonding combinations of transition metal d and 
oxygen p states. In the single crystal form of many tran- 
sition metal oxide materials the transition metal-oxygen 
hybridization vanishes at the zone center (F) point of 
the Brillouin zone. The bands are degenerate or nearly 
degenerate at this point, and disperse upwards from it. 
A partial breaking of translational symmetry activates 
the hopping at zone center, thus lifting some orbitals up 
relative to others. Chaloupka and Khalliulin^ recently 
proposed that in a heterostructure composed of alternat- 
ing layers of LaNiO^ and LaAlO^ this effect could lead 
to a situation in which only the — Ni Cg symme- 
try orbital would be occupied, providing an electronic 
structure similar to that found in high Tc copper-oxide 
superconductors. The issue was studied theoretically by 
Hansmann and coUaboratorsiS. 

In this paper we identify a third route to control of 
orbital polarization, based on chemical composition of 
spacer layers in an oxide superlattice. We consider specif- 
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ically the class of systems introduced by Chaloupka and 
KhaliuUin^, namely superlattices composed of alternat- 
ing layers of LaNiOs and of a spacer layer LaXOs the or- 
bital polarization depends upon the choice of counterion 
X, even when X is chosen such that LaXO^ is a wide- 
bandgap insulator. In this situation the conventional wis- 
dom is that because the LaXO^ layers are insulating, the 
only relevant effects are strains induced by lattice mis- 
match between the LaNiOs and the LaXOs layers. We 
show that this is not the case, and that the strength of 
chemical bonding of the apical oxygen with the X ion sig- 
nificantly affects the polarization, even changing its sign 
in some cases. 

The rest of this paper is organized as follows. Sec- 
tion [IT] presents the systems wc study, the definitions we 
use and the formalism we apply. Section IIIII presents 
the band structures and resulting orbital polarizations in 
the simple case where lattice relaxations are forbidden. 
Section IIVI presents a tight binding analysis which expli- 
cates the physics behind the results presented in Section 
IIIII Section |V] considers the additional effects of lattice 
relaxation. Section IVII is a summary and conclusion. 

II. MODEL, DEFINITIONS AND METHODS 

We consider superlattices of the form depicted in Fig. 
[1] composed of blocks of LaNiOs and another material 
LaXO^, alternating along the (001) axis of the basic 
perovskite crystal structure. Wc choose X=B, Al, Ga, 
In from column 3A of the periodic table. While to our 
knowledge only LaAlO^ has been used to grow oxide su- 
perlattices, all of the LaXOs compounds have been re- 
ported in the literature^ii^iii. For simplicity we take all 
Ni-O-Ni and Ni-O-X-O-Ni bonds angles to be 
180°. In most of our calculations we fix the in-plane lat- 
tice constants to the GGA-optimizcd pseudocubic bulk 
LaNiOs value 3.8lA but in some of our calculations the 
in-plane lattice constant was set to 3.9lA to model the 
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FIG. 1: Central portion: Unit cell of (001) superlattice con- 
sidered in this paper with atoms indicated on right. Left 
side: chain of hybridizing orbitals (top to bottom Xs, Op. 
and Ni^^i_^i) controlling relative occupancy of 0(322 or- 
bital. 
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FIG. 2: d-projected density of states (upper panel — 
lower panel — '■'^) of (001) La2NiXO(i heterostructures 
with counter-layer B-site ion X—B (black) , Al (red), Ga 
(green) and In (blue), obtained from GGA band theory cal- 
culations for unrelaxed structure (ideal cubic perovskite, dou- 
bled in (001) direction) with all Ni - O and X - O bond 
lengths taken to be one half of the bulk LaNiOz lattice pa- 
rameter a — 3.8IA. Arrow in lower panel indicates region 
where dependence of d^^2_j.2 density of states on counterion 
is evident. 



effect of a film grown on an SrTiO^ substrate. We con- 
sider two cases, an 'unrelaxed' case where the out of plane 
lattice constant is set equal to the in-plane one, and a 're- 
laxed' structure where the out of plane atomic positions 
are adjusted to minimize the energy as described below. 

To determine the electronic structure we use density 
functional band theory within the Generalized Gradi- 
ent Approximation (GGA) as implemented in the Vi- 
enna ab-initio simulation package (VASP)ii"— . Wc used 
a plane- wave basis set and the projector augmented wave 
methodic with a cutoff of 270 eV and fc-point meshes of 
10 X 10 X 5. In the 'relaxed' calculations, atomic positions 
were relaxed along the (001) using conjugate gradient 
minimization of the GGA energy. 

To determine orbital occupancies we projected the cal- 
culated electronic density of states onto locally defined 
atomic orbitals obtained by defining a sphere around the 
atom in question and then projecting the wave functions 
within the sphere onto the appropriate symmetry states. 
The sphere sizes in A were taken to be 1.286 {Ni-d) 0.820 
{O-p) (these are the VASP defaults) while for the X-s or- 
bital on all of the countcrions wc used 1.402. 

Examples of our results for d-projected densities of 
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FIG. 3: Integral of Ni-do,^i_^i (red solid line) and Ni-d^2_y2 
(blue dashed line) density of states for the four choices of X 
ion considered in the paper. Flat region between ~ —2 and 
~ — 1 eV defines lower edge of conduction band. 



states are shown in Fig. [2l The density of states consists 
of two components: a broad antibonding band (of mixed 
Ni-d/O-p character) spanning the region from —1.5 to 
^ 2eV near the chemical potential and a narrow bond- 
ing band at the low energy E —6eV. As can clearly be 
seen, the lower edge of the Ni — O antibonding band is 
well defined. Orbital occupancies are then obtained from 
integrals of the densities of states; representative results 
are shown in Fig. |3l That the integrals are not quite 
2 is related to the sphere size. We are interested in the 
occupancies nx2_y2 and n^z^-r'^ of the near fermi level 
d — p antibonding bands. To obtain these we take the 
difference between the value of the density of states inte- 
gral at the lower edge of the antibonding band (identified 
from the flat part of the integrated density of states plot) 
and the value at the fermi level. Values corresponding to 
integrating over the entire conduction band plots of inte- 
grated d spectral weight vs energy may be read directly 
from Fig [3] and yield the same conclusions. 

From the orbital occupancies we obtain the orbital po- 
larization P which we define as 



III. RESULTS: UNRELAXED STRUCTURE 

Table [J presents our computed results for the orbital 
polarization. In this section we focus on the second col- 
umn, giving results for 'unrelaxed' structures in which 
all Ni — O and X — O bonds set equal to one half of 
the Ni — Ni distance of bulk LaNiO^. Results from 
these structures highlight the chemical effect of interest 
here. Results for 'z-relaxed' structures obtained by min- 
imizing the energy with respect to atomic motions along 
the (001) (supciiattice) direction and for structures with 
additional in-plane strain are discussed in section Ivl 
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.26 
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TABLE I: Orbital polarization calculated as defined in Eq ^ 
for three superlattice families: 'unrelaxed LNO', with all A'^i — 
O and X — O bond-lengths set equal to the GGA-optimized 
pseudocubic LaNiOz value of 1.905A ; 'z-relaxed LNO', with 
in-plane Ni — O and X — O bonds set equal to 1.905A and 
out of plane bonds relaxed to minimize the GGA band theory 
energy and 'relaxed STO' in which the in-plane Ni — O and 
X — O bonds are set equal to 1.95 A and out of plane bond 
lengths are relaxed to minimize the GGA band theory energy. 

The calculated polarizations are seen to vary strongly 
with choice of counterion X. The polarization differences 
can also be seen directly in Fig. [21 for example as an X- 
dependent change in the d^^2_j.2 density of states in the 
lower portion of the antibonding band (indicated by ar- 
row). We stress that the polarization differences occur 
even though all Ni — O and X ~ O bond lengths are 
equal; thus the difference is a chemical effect. We also 
stress that while the magnitude of the polarization de- 
pends on calculational details such as the sphere sizes 
and the range over which one integrates, the trends be- 
tween materials are robust and clearly demonstrate that 
a non-structural difference between the different X ions 
strongly influences the orbital polarization. 

We believe that this non-structural difference is re- 
lated to the properties of a near-fermi-level orbital on 
the X site. To investigate this hypothesis we computed 
orbitally projected densities of states for the different sys- 
tems. We found the most significant effect comes from 
the s-symmetry orbital on the counterlayer X site. In 
Fig. |4] we present the s projected density of states on 
the X site, along with the Ni — and apical oxy- 
gen pz projected densities of states. Examination of the 
series B, Ga, In indicates a clear correlation between 
the orbital polarization and the energy of the X-site s- 
symmetry orbital and its hybridization with the apical 
oxygen orbital. The Al case is an outlier in this series, 
for reasons not yet understood. 
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FIG. 4: Orbital symmetry-projected densities of states for un- 
relaxed La2NiXOe heterostructures with X=B (top panel), 
Al (second from top), Ga (third from top) and In (bottom). 
Shown are the — orbital on the Ni (green trace), the 
Pz orbital on the apical oxygen (blue trace) and the s orbital 
on the X site (red trace). 



ultimately to changes in the hybridization between the 
transition metal ion and the nearby oxygen orbitals. In 
the La2NiXOe, systems the crucial role is played by the 
apical oxygen which as can be seen from Fig. [T] connects 
the s and Ni and couples to the Ni — orbital but 
not the — orbital. Fig. U shows that the different 
choices for ion X have states with differing overlap with 
the Pz state on the apical oxygen. Thus in essence the 
counterion affects the polarization by shifting the prop- 
erties of the apical oxygen pz state. 



IV. TIGHT BINDING ANALYSIS 

To confirm the hypothesis that the trends reported in 
Fig. |4] are causally related to the observed polarization 
changes we turn to a tight binding model which, while 
necessarily a simplified description of the actual band 
structure, captures with reasonable fidelity the essential 
features of the density of states. The specifics of the tight 
binding model are given in the Appendix. We present 
here the main ideas. We begin with a 5-band model of 
cubic LaNiO^ consisting of the two Cg symmetry Ni d- 
orbitals and the three oxygen p^ orbitals. We find that 
a very good representation of the GGA-calculated Ni eg 
and Opa density of states for cubic LaNiOs is obtained 
with a Ni d- level energy of —l.22eV, an oxygen energy 
of —5.2eV, and Ni — O and O — O hopping amplitudes 
of l.SeV and 0.7eV respectively. We next double the 



Higher orbital polarization {In case) is associated with 
an s orbital which is farther removed in energy and less 
strongly admixed with Opz- Higher orbital polarization 
is also associated with stronger mixing between the Ni 
3^2 _ j,2 orbital and the apical Op^ . Lower orbital polar- 
ization (B case) is associated with an s orbital which is 
closer to the fermi level in energy and more strongly ad- 
mixed with Op^ and is also associated with weaker mixing 
between the d^z^-r^ and the apical O. 

Our finding of a key role played by the s orbital on the 
X site is reminiscent of results of O. K. Andersen and 
collaborators These authors argued that in the high 
Tc cuprate case the variation of the fermiology across 
different sub-families of cuprate materials was controlled 
by the energy of the Cu 4s orbital, which affected the 
ratio of first and second nearest neighbor hopping t' /t 
(similar arguments relating the fermiology of orbitally 
polarized nickelate heterostructures to the Ni — 
orbital were made by Chaloupka and Khalliulin^) . Here 
we argue in an analogous way that the variation of the 
energy of the X-site s orbital controls the polarization 
of the Ni d-orbitals. (We also investigated the Ni As 
orbital, finding that its energy does not vary significantly 
across the series). 

We suggest that the polarization differences are due 
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FIG. 5: Density of states computed from tight binding model 
projected onto Ni d3,2_^2, apical Op^ and Xs states. 
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unit cell in the z direction and replace one of the two Ni 
with an X orbital. We must then introduce three new 
tight binding parameters: an energy ex of the orbital on 
the X site and hopping parameters tsp,tspz coupling 
the X orbital to the oxygen ions in the X plane and to 
the apical oxygens. The parameters ex, tsp, and tspz 
are numerically optimized by minimizing the difference 
between the GGA and tight binding (TB) densities of 
states of Ni^^2_j.'2^ apical Opz and Xs in the near fermi 
level and positive energy regions. Results are shown in 
Fig. [5l We see that the tight binding model reproduces 
the behavior of the near fermi level states well. Differ- 
ences (for example an overestimate of the width of the 
low energy E ^ —12eV portion of the s band and an 
underestimate of the width of the high energy E ^ 8eV 
portion) arc present, but total number of X-s states in 
the lower and upper energy portions agree well with the 
GGA calculations (see Appendix). 

The orbital polarizations were then computed from the 
tight binding model, with the results Pp = 0.19, Pai = 
0.28, Pca = 0.31 and P/„ = 0.39. Comparison to the 
information presented in Table [J shows that the tight 
binding polarizations reproduce well the trends found in 
the GGA calculation. Because the only 'moving part' in 
the tight binding model is the X orbital the successful 
fitting confirms that the variation in orbital polarization 
between materials is in fact due to changes in nominally 
unfilled orbitals on the X site. 



V. EFFECT OF STRUCTURAL RELAXATIONS 

We finally consider the interplay of structural relax- 
ations with the effects we have studied so far in this 
paper. We performed a structural relaxation process in 
which atoms were allowed to move in the direction trans- 
verse to the plane so as to minimize the GGA energy 
with the in-plane lattice constants kept equal to the bulk 
LaNiOs value. Table |TT] presents the TVi-apical O and 
X-apical O bond lengths found after relaxation. Note 
that the TVi-apical O bond length is remarkably robust: 
almost all of the bond length change due to relaxation 
occurs in the X-O bond. We take this as further confir- 
mation of the 'chemical effect': the geometrical Ni — O 
distance is less important than what happens at the other 
end of the Ni^O~X bond. The corresponding polariza- 
tions are given in the third column of Table IH Relaxation 
has a large effect on P, and acts to enhance the 'chemi- 
cal' effects we have identified. Note in particular that the 
sign of the polarization actually reverses in the B case. 

The first and third panels of Fig [6] show the GGA den- 
sities of states of the two extreme cases {B and In) com- 
puted for the relaxed structure. Comparison to Fig. [3] 
shows that in the boron case the Xs density of states 
shifts upwards in energy, while in the In case the Xs 
density of states shifts downwards in energy, as expected 
if the position of the X orbital is controlled by hybridiza- 
tion with the oxygen pz . 



X 


jNi-O ,X-0 
"■LNO "-LNO 


,Ni-0 ,X-0 
"STO "STO 


B 
Al 
Ga 
In 


1.07 0.85 
1.00 1.00 
1.00 1.06 
0.99 1.19 


1.05 0.80 
0.99 0.98 
0.99 1.05 
0.97 1.16 



TABLE II: Bond lengths between Ni and apical oxygen and 
between X-site atom and apical oxygen presented as ratio of 
computed bond length to LaNiOa value a = 3.81A for relaxed 
structures with in-plane lattice constant set to LaNiOs value 
(LNO) and to SrTiOs value a = 3.905 (STO). 
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FIG. 6: GGA orbital symmetry-projected densities of states 
for relaxed La2NiXOe heterostructures with X=B (top 
panel), and In (third from top), along with tight binding 
fits for the two cases (second and fourth panels from the top). 
Shown are the Sz^ — orbital on the Ni (green trace), the 
orbital on the apical oxygen (blue trace) and the s orbital 
on the X site (red trace). 



To show how the effects of structural relaxations are 
incorporated within the tight binding model we have used 
the same TB parameters as in the unrclaxcd case, except 
that we have increased tspz in the B case (modelling 
the effect of a decreased X — O distance and decreased 
it in the In case (modelling the effect of an increased 
X — O distance. The resulting tight binding densities 
of states are also shown in Fig. [6] The tight binding 
model reproduces the basic shifts in the X orbital density 
of states and leads to polarizations, Pg^ = 0.08 and 
Pf„«= 0.44. 
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Finally, to study the effects of in-plane strain (induced, 
for example by growing on a different substrate) we set 
the in-plane lattice constant equal to the value 3.9lA ap- 
propriate to the widely-used SrTiO^ substrate material. 
This corresponds to applying a tensile in-plane strain. 
Results for the polarization are shown in the fourth col- 
umn of Table H] and for the lattice constants in the fourth 
and fifth columns of Table [TTl The tensile strain affects 
the polarization in the expected way but the out of plane 
lattice constant is again remarkably robust. 

VI. CONCLUSION 

Thus in summary we have shown that in addition 
to the well-established geometrical effect, orbital polar- 
ization in an oxide heterostructure may be controlled 
by appropriate choice of counterions in the superlattice. 
Changes in the counter-ion produced large orbital po- 
larization differences, which are expected to be further 
enhanced by correlation effects (not yet included in our 
calculations). Both the magnitude and (after relaxation) 
the sign of the polarization can be changed. The crucial 
factor was found to be the hybridization of an s-like or- 
bital on the counter-ion with the Pz orbital on the apical 
oxygen. This hybridization changed the interaction of 
the apical oxygen pz with the transition metal ion, and 
hence changed the orbital polarization. We also found 
that under structural relaxation the Ni — O bond lengths 
change relatively little, while the X — O bonds changed 
substantially. Our results thus show that the transition- 
mctal/oxygen bond length is not the only variable con- 
trolling polarization, and that chemical effects should in 
general be considered when attempting to optimize su- 
perlattice properties. We remark that the need to in- 
clude the s orbital demonstrates a limitation of the com- 
mon theoretical strategy of deriving from band structure 
a 'minimal low energy model' of the system of interest. 
While of course a careful downfolding, keeping track of 
the effects stemming from bands that are projected out, 
will produce a low energy model which contains the ef- 
fects of interest (see, e.gi "*-^!^^ ), construction of 'minimal 



Here I3 is the 3x3 unit matrix and we took ep = 
—5.2eV and tpp = 0.7eV. Inclusion of oxygen-oxygen 
hopping is necessary to reproduce the narrow d feature 
observed at ~ —QeV in e.g. Fig. 2 of the main text. 
Ho is just Hp doubled in the z direction. 



high energy models' such as the one we have defined may 
also be a useful strategy. 
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Appendix: Details of tight binding model 

The tight binding model presented in the text is based 
on the minimal assumption of a cubic ABO3 perovskite 
structure, with the unit cell doubled in the (001) direction 
and the two B-sites distinguished. It involves 9 orbitals: 
Ni — y'^ and 3z^ — r'^, the Op^ orbitals in the Ni-0 
plane {Oni.x, ONi,y), Opa orbitals [Oapi.z, Oap2.,z) at the 
two apical oxygen sites, an s orbital and two Op^ orbitals 
in the X — O plane {Ox,x, Ox,y)- 

Ordering the basis states as (3z^ — r^,x^ — 
y'^,ONt,x,ONt,y,OAPi,z,Ox,x,Ox,y,OAP2,z,X) wc may 
write the tight binding Hamiltonian Htb in a schematic 
block-diagonal form as 

/ Hm hI^ \ 
Htb = HpD Ho H^^ (2) 
\ Hpx Hx J 

We took Hj^i = ef/l2with / the 2x2 unit matrix and 
Ed = -1.22eV. 

For the oxygen portion Ho we assumed the usual over- 
laps between orbitals on second neighbor oxygens, wrote 
the form appropriate for a simple cubic lattice, and then 
doubled the unit cell along (001). To avoid a cumber- 
some display of 6 x 6 matrices we write the result before 
doubling as 



- e-'*-") (1 - e'''-) (1 - e"**^-) \ 

(1 - e'*^") (1 - e-*'^-) (3) 
-e-**^") J 

I 



For the p — d hopping we took the standard form 




Hr, 



1 - e-''^- 

-ife. 



— e 



1^1 - e^'^y] 



(1 



(1 
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with tpD = l.SeV^ while for Hxp we put 



Hpx — 






tsp (l - e 



\ 



(5) 



The important parameter is tspz giving the overlap 
between the s orbital and the apical oxygen states. 

The tight binding parameters are not uniquely deter- 
mined; changes in one parameter may be to some ex- 
tent compensated by changes in another, but as long as 
the density of states and s-apical Pz mixing arc repro- 
duced with reasonable accuracy the polarization is ro- 
bust. We chose to fix the tight binding model param- 
eters by first choosing the iVi-related parameters to fit 
the GGA band structure of cubic LaNiO^. Then we 
determined the other orbitals by minimizing the 'dis- 
tance' between the GGA and TB predictions for the 
density of states of the Sz^ — r^, Opz and Xs orbitals. 
That is, we chose the TB parameter set which minimizes 
J2o I ^^iP(^)GGA ~ p{^)tb)'^ where p{e)'-' denotes the 
density of states for orbital O = Ni d^z^-r'^, Opz or 
Xs atom and the energy ranges arc —1 to 2 for the Ni, 
— 10 to the upper-limit for the Opz and —2 to the upper 
limit Xs orbitals. The parameters used to produce the 
fits shown in the text for the unrelaxed structures are 



summarized in Table IIIII along with the resulting polar- 
izations. To model the effects of structural relaxation in 
the In case we decreased the tspz from 4 to 2.5 and in 
the B case we increased it from 4.5 to 6. The resulting 
polarizations are Pj^ = 0.08 and P^^ = 0.44. 





tsp 


tspz 


ex 


P 


B 


3.0 


4.5 


-4.0 


0.19 


Al 


5.0 


6.0 


-1.6 


0.28 


Ga 


4.5 


5.0 


-1.6 


0.31 


In 


4.5 


4.0 


3.2 


0.39 



TABLE III: Tight binding parameters used to construct fits 
displayed in Figs. 4 of main text, with resulting orbital po- 
larization 



Figure [7] compares the integrated X-s density of states 
obtained from GGA and optimized TB calculations re- 
spectively. Some differences are visible. At the very low 
end of the energy range { E ^ io —12eV the X den- 
sity of states extends to low energies in the tight binding 
model than it does in the GGA calculation, presumably 
because of level repulsion arising from other low energy 
orbitals present in the GGA calculations but not in the 
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FIG. 7: Left panel: integral of X-s GGA density of states for the four choices of X ion considered in the paper. Right panel: 
integral of X-s TB density of states for the four choices of X ion considered in the paper 



